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We show how primordial bispectrum of density perturbations from inflation may be characterised 
in terms of manifestly gauge-invariant cosmological perturbations at second order. The primordial 
metric perturbation, f, describing the perturbed expansion of uniform-density hypersurfaces on 
large scales is related to scalar field perturbations on unperturbed (spatially-flat) hypersurfaces at 
first- and second-order. The bispectrum of the metric perturbation is thus composed of (i) a local 
contribution due to the second-order gauge-transformation, and (ii) the instrinsic bispectrum of 
the field perturbations on spatially flat hypersurfaces. We generalise previous results to allow for 
scale-dependence of the scalar field power spectra and correlations that can develop between fields 
on super-Hubble scales. 
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Q 1 I. INTRODUCTION 

<d : 

Inhomogeneities in the distribution of matter and radiation in the universe today provide valuable information about 
■ the dynamical history and physical processes in the very early Universe. In particular temperature anisotropies in the 
t-H | cosmic microwave background (CMB) sky provide detailed information about the primordial density perturbations 
at the time of last scattering of the CMB photons. At present the simplest explanation for the origin of these 
perturbations seems to be that they originated as quantum fluctuations in one or more light fields during a period 
of inflation in the very early Universe pi Q] • These vacuum fluctuations on small (sub-Hubble) scales were stretched 
to large (super-Hubble) scales by the accelerated expansion, where they subsequently evolved as effectively classical 
\ perturbations. 

Q\ ■ Vacuum fluctuations of massless scalar fields in a de Sitter geometry enter a squeezed state on scales larger than 
the Hubble length where the decaying mode can be neglected. Thus they can be treated as an effectively classical 
. distribution on large scales with Gaussian statistics. For a Gaussian random field the power spectrum is sufficient to 
describe all the statistical properties of the distribution. But non-linearities in the evolution of the initial fluctuations 
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will lead to a non-Gaussian distribution. Until recently studies of non-Gaussian perturbations from inflation were 
largely restricted to self-interactions of scalar fields |j, |a| neglecting non-linear gravity. Gravitational effects can be 
included in a stochastic approach using local equations to study the evolution of the coarse-grained inflaton field on 
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2 1 long wavelengths [fil LI |8|, |9| . 

Non-Gaussianities inevitably arise at second-order in cosmological perturbation theory but the nature of the non- 
Gaussianity will depend upon the choice of variables [T(i| . Acquaviva et al constructed a gauge-invariant quantity 
at second-order on large scales , bu t its physical interpretation was left open and it was only approximately constant 
during slow-roll inflation 0, Maldacena gave an analysis of second-order perturbations from slow-roll 

1 inflation working in terms of the field perturbation on unperturbed spatial hypersurfaces and relating these to the 
curvature perturbations on uniform field hypersurfaces during single field inflation. Both these quantities have a clear 
physical definition and hence are implicitly gauge-invariant. He showed that there is a curvature perturbation, which 
becomes constant on large scales and whose bispectrum for squeezed triangles (k% <C ^2,^3) is proportional to the 
tilt of the scalar power spectrum, and hence small. Maldacena's result has subsequently been verified by a number of 
authors [Til flEl iwL ITtI fl8| , going beyond slow-roll , and extending it to loop corrections an d non- minimally 
coupled fields |15^ For a recent review see Ref. [2(j ■ 

Onl y m ultiple field models of inflation seem likely to generate significant non-Gaussianity in the density perturbation 
|M ll^ l2l(. 123. 123 . 12^ . l25l l26j . In particular Lyth and Rodriguez |22| have recently emphasized that the non-Gaussianity 
of the curvature perturbation may be simply written down in the so-called "JTV-formalism" [23, l28l I2H |3(j where 
one uses the "separate universes" approach |31| to follow the non-linear evolution on super-Hubble scales in terms of 
locally homogeneous solutions. 

In this paper we show how the non-linear evolution of primordial perturbations from inflation may be characterised 
in terms of manifestly gauge-invariant perturbations at second-order. We adopt the formalism of Malik and Wands 
|32j who showed that there is a conserved curvature perturbation at second-order for adiabatic density perturbations 
on long wavelengths. In section II we define the gauge-invariant curvature perturbation at second-order and give 
the leading order expressions for the power spectrum and bispectrum. In section III we specialise to the case of 
adiabatic perturbations on large scales from single-field inflation. We relate the conserved curvature perturbation 
on large scales to the gauge-invariant field perturbations on spatially flat hypersurfaces during slow-roll inflation, 
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consistent with the <5iV-formalism. In particular we distinguish the purely local form of the bispectrum generated 
by the second-order gauge transformation from the intrinsic bispectrum of the field perturbations on unperturbed 
spatial hypersurfaces. We use a simple argument to estimate the bispectrum of field due to large-scale variations in 
the local Hubble rate during inflation, and hence recover Maldacena's result for the non-linearity of the curvature 
perturbation. We generalise these results to the multiple-field case in Section IV, allowing for scale dependence of 
the power spectra and cross-correlations between the fields which may develop on large scales. These results can be 
simplified by decomposing the field perturbations into an infiaton component (along the instantaneous background 
trajectory) and orthogonal isocurvature perturbations [33J. We conclude in Section V. 



II. GAUGE-INVARIANT PERTURBATIONS 



Discussions of relativistic perturbations have historically been plagued by gauge-dependence of apparently physical 
quantities such as the density perturbation. For instance, in a spatially homogeneous background where the density 
evolves in time, the inhomogeneous density perturbation depends on the choice of time coordinate at first-order, and 
spatial coordinates at second-order. This "gauge-problem" can be avoided by using a physical choice of gauge and 
building gauge-invariant definitions of the perturbations in that gauge. Thus Bardeen constructed gauge-invariant 
cosmological perturbations at first-order [24j and Malik and Wands have shown how to construct gauge-invariant 
cosmological perturbations at second-order 32]. This still leaves a number of different possible gauge- invariant 
definitions of the density perturbation, for example, corresponding to different physical gauge choices, though the 
perturbations are nonetheless gauge- invariant. 

We will split any scalar (f> into a spatially homogeneous background and first- and second-order inhomogeneities: 

cj)(t, x) = <p (t) + Si^>(t, x) + -5 2 <p(t, x) . (1) 

The split between first and second order inhomogeneities for any given initial perturbation is arbitrary, so we can 
use this freedom to specify that the first-order perturbations are Gaussian random fields and any non-Gaussianity 
is described by the second-order perturbation. In particular we will take the Gaussian part of the perturbations to 
originate from free field fluctuations at early times during inflation, and assume the second-order perturbations are 
vanishing in this early time limit. 

Under a temporal gauge transformation, a = a± + (0:2/2), a scalar (j> transforms at first- and second-order as 

5i4> — > <5i0 + 0oai , (2) 
S 2 (f> — * S 2 cf> + 0o"2 + d t ^oai + 2(5i^ a\ . (3) 

We will adopt the method of Malik and Wands [s3 to define gauge-invariant cosmological perturbations up to second 
order by constructing a physical perturbation in a physically defined coordinate system. 

There are (at least) two different ways of writing the perturbed line element on any spatial hypersurfaces: 

ds 2 = a 2 (i) [(1 - 2i))8 lJ dx l dx j + 2E tij + 2F (iJ) + fty] , (4) 
= a 2 {i) [e 2SN 5 ij dx i dx j + + 2F (U) + h i3 ] , (5) 

where Fi is the divergence-free, vector perturbation, and hij is the transverse and trace-free tensor perturbation. One 
may refer to ip as the curvature perturbation and SN as the perturbed expansion, though the two are trivially related 
as -0 = (e 2SN — l)/2. To first and second order we write 

tp = i>i + . (6) 

SN = S-^N+^N, (7) 

and hence 

V>i = -SiN, (8) 

02 = -S 2 N - 2(5 1 N) 2 . (9) 

The curvature perturbation, ip or SN, on uniform-density hypersurfaces is non-linearly conserved on large scales 
(where gra dient terms can be neg lected) for adiabatic perturbations as a direct consequence of local energy conservation 
|29l I30I l3ll 13^] . (see also [35:, 36]). Thus we will be particularly interested in calculating this perturbation. 
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Most authors [l3L 12^. I301 137^ have worked in terms of SN on uniform-density hypersurfaces, which we denote by 



Ci = 5 1 N-—S 1 p, 
Po 

J- — X AT H , J_ ( H5 ^P 

Q 2 = o 2 N-—b 2 p+\ — : — 
Po \ Po 




(10) 
(11) 



This definition leaves a residual dependence at second order on the choice of spatial gauge, i.e., the choice of spatial 
coordinates threading the uniform-density hypersurfaces. This can be eliminated by a physical choice of spatial gauge 
at first-order |32l |. such as working with worldlines orthogonal to the uniform-density hypersurfaces. However in what 
follows we will consider perturbations on large scales where in any case the residual spatial gauge dependence becomes 
negligible 29]. 

Malik and Wands [32^ gave a gauge-invariant definition of the curvature perturbation at second order on uniform- 
density hypersurfaces, ip\ p , and showed this is conserved for adiabatic perturbations on large scales. This is identical 
to — £ defined above at first-order, tpx\ p = — £i, but differs at second order simply due to the difference in choice of 
curvature or expansion perturbation @. Thus we have |10| 



(12) 



Hence both £ and ip\ p are conserved at second order for adiabatic perturbations on large scales. 

For a stochastic random field, as would be expected to arise from vacuum fluctuatuations during inflation in the 
early universe, the perturbations are most usefully described in terms of their power spectrum, and higher-order 
moments, in Fourier space. At linear order the Fourier coefficients of ( are defined by the same equation (|10|) as in 
real space, while at second order we have 



C2.k = S 2 N k - —S 2Pk + — 
Po \ V Po 




where "*" represents a convolution 



/DO 
/(k-k')<?(k') 
-oo 



(13) 



(14) 



where k represents d 3 k/(2ir) 3 . We will be using the corresponding notation (F 3 (k) to represent (27r) 3 5 3 (k) 
The power spectrum of the perturbation is given by 

p c (k) ^(k + kOMCkCO, 



(15) 



For an isotropic distribution the power spectrum is a function solely of A; = |k|. Note that the variance per logarithmic 
interval in fc-space is given by 



V c (k) 



4nk 3 
(2^F 



P C (k). 



The bispectrum of the curvature perturbation is 

B c (ki,k 2 ) ^ 3 (k x + k 2 + k 3 ) = (Ck x Ck 2 Ck 3 ) , 

which by the requirement of statistical isotropy is non-zero only for a triangle of modes, ki 
The bispectrum is zero for Gaussian perturbations and hence to leading order we have 

p c (k) J 3 (k + k0 = (Ci,kCi,0, 

B C (ki»k a ) ^ 3 (k 1+ k 2 + k 3 ) = i[(Ci,kCi,k'C 2 ,k") 

+ (Cl,kC2,k' Cl,k") + (C2,kCl,k' Cl,k")] 



k 3 = 0. 



(16) 



(17) 



(18) 



(19) 



If the second order part of the curvature perturbation can be written purely in terms of the square of the first order 
part in real space 



C2 



-(6/W5)Ci 



(20) 
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where /at l is a function of background parameters then the bispectrum will be given by 

B c {k u k 2l k 3 ) = -^fNL(Pc(ki)P C (k 2 ) + P f (fci)P C (fc3) + Pdk 2 )Pc(k 3 )) ■ 



(21) 



calculated for closed triangles, ki + k.2 + k3 = 0. Jnl is then a quantifier of this form of "local" non-Gaussianity and 
is called the non-linearity parameter. 

III. SINGLE-FIELD INFLATION 

We will now calculate the bispectrum for the scalar field perturbations on large-scales in single-field models of 
inflation, and hence reproduce the bispectrum for the large-scale curvature perturbation, £, [l3j . 

A. Curvature perturbations from inflaton perturbations 

During single field inflation, driven by field tp, density perturbations become adiabatic in the large-scale limit at 
first- |33| or second-order . An important consequence is that ( is conserved and the uniform-density hypersurfaces 
coincide with uniform-field hypersurfaces on large scales. Thus we have 



& 

C2 



8 l N-—5 l ip, 
<P 

S 2 N - ^-S 2 ip 

tp 



tp 



Sup SxN 



P 



H 



Hi 




Sip 




-i] 





(22) 
(23) 



These expansion perturbations on uniform-field hypersurfaces are closely related to the field perturbations on 
uniform-expansion hypersurfaces, which have the gauge-invariant definitions 



5 X Q = 5i<p-?=8iN, 

ti 

x n - x V x at , ( <P 6lN 
S 2 Q = 5 2 p - —S 2 N + I — 



<5iA Si(p 



H 



tp 



^~h) \h~ 



(24) 
(25) 



At first-order S\Q is the usual Mukhanov-Sasaki variable |38l l39| used to q uantise the semi-classical linear perturba- 
tions during inflation, and 5 2 Q is its natural extension to second-order |13ll32|. 

We can re-write the metric perturbation £ in Eqs. (|22ll and (|23|l somewhat more succinctly in terms of these 
gauge- invariant field perturbations, to give (see also Ref . |40j ) 



Ci 

C2 



H, 

<P 



-SiQ. 



-^rS 2 Q 

tp 









\h 2 


Hip 





tp 



(26) 
(27) 



where we have used 



SjN Sxtp = Ci 
H tp H 



H 



H 2 Hip 



<P \ Slip 



tp 



(28) 



Introducing the slow-roll parameters 



H 

\H 2 Hp> 



(29) 
(30) 
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and using the fact that (i = for adiabatic perturbations on large scales, we finally have 

G = --S 2 Q + (2e- V )(^-) 2 . (31) 

<P \ 9 J 

We can show that the non-linear adiabatic curvature perturbation on large scales during inflation can simply be 
written as the perturbed expansion history, N = J Hdt, due to the gauge-invariant field perturbations [2^, [2!| 



This is a non-linear extension of the <57V-formalism 27, 28] for calculating perturbations from inflation. Hence to first 
and second order we have 

Ci = f^iQ, (33) 



where 

dN H 

dip <p 



(35) 



-^r r- =(2e- V )r- . (36) 



d 2 7V ( H (p\ (H\ 2 (H^ 2 



dip 2 \ H 2 Hip I \p J \tp 



B. Non-Gaussianity of field perturbations 

During slow-roll inflation the initial state of the scalar field perturbations can be set by imposing the Minkowski 
vacuum state on small scales (much smaller than the cosmological Hubble scale, k/a 3> H) for the scalar field 
perturbations on spatially-flat hypersurfaces, SN — 0. For any weakly-coupled scalar field this gives an effectively 
Gaussian distribution of initial field fluctuations. For any light field (with effective mass less than the Hubble scale) 
living in a de Sitter spacetime, the Hubble damping leads to a squeezed state on scales larger than the Hubble length 
and hence we can treat the field as a classical random field with 

^2(M = !J, (37) 

and hence V(;(k*) = (H/2tt) 2 , where a * subscript denotes Hubble-crossing, fc* = aH. 

We will give a s imple extension of the argument due to Maldacena 14] to calculate the amplitude of the second- 
order (non-Gaussian) contribution to the scalar field perturbation due to gravitational coupling at Hubble-exit. 

It will be convenient to write the first-order field perturbation at Hubble exit as 

StQ(K) = ^e k , . (38) 

where e k is a classical Gaussian random variable with unit variance: 

/A „ . (F(k-fk') 

{eiAf ) = t±^2. (39) 

However at second order there is a local correction to the amplitude of vacuum fluctuations at Hubble exit due to 
first-order perturbations in the local Hubble rate, H. This is determined by the local scalar field value due to longer 
wavelength modes that have already left the horizon 

H = H(f Q ) + H\p Q ) [ C d 3 k 6iQ k . (40) 
Jo 



where k c < fe* is a cut-off wave-number. 
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Thus for fc* 3> k c we have upto second order 

<>Q(k*) = ^-e k + / fk^^e k ^ k . (41) 

Z7T J Z7T 

and hence we can identify the second order scalar field perturbation at Hubble exit: 

U 2 Qy* = ^(<5i<9»5iQ)k* 5 (42) 
Z ti 

where we have taken <5iQk„-k' = (-ff/27r)e k ,- k ' f° r ^* 3> an d a "•" represents a cut-off convolution, that is a 
convolution-type integral defined by 

(f»g) k = f ° d 3 k'f k ,g k _ k , , (43) 
Jo 

where fe c -C fc. 

This effect gives a non-zero bispectrum ()17|l for the scalar field perturbations. In the squeezed triangle limit, k 2 — k^ 
and k\ <C k 2 , the correlation of the field perturbations at the point when the two smaller modes cross the horizon 
(&2 = o,H) gives, to leading order, 

(6Q kl 5Q k2 SQ ks ) = {S!Q kl ^5 2 Q k2 <SiQk 3 ) + (5iQ kl <5iQ k2 ^2<3k 3 ) , (44) 

Note that the second order part of the perturbation on the larger scale, S 2 Q kl , will be uncorrelated with the modes 
on smaller scales. Substituting our expression (|42[1 for the second-order field perturbation at horizon crossing into 
Eq. yields 

(SQ kl SQ k2 5Q k3 ) = 2^(5 1 Q kl 6 1 Q k2 (6 1 Q»5 1 Q) k3 ), (45) 
= -2e (j) P Q (A; 1 )F Q (fc 2 ) ^ 3 (k x + k 2 + k 3 ) (46) 
Comparing this with the definition of the bispectrum in Eq. (|17f) we see that 

B Q (k,h)=-2e(^]p Q (k)P Q (K). (47) 

Note that this is not strictly of the local form given in Eq. I|21|) due to the asymmetry between long wavelengths, 
k\, and shorter wavelengths, k 2 and k%. However P(k) P(k*) for k <C fc* as V((k) oc k 3 P^(k) generated during 
slow- roll inflation is nearly scale- invariant, and hence this difference is not significant in the squeezed triangle limit. 

One can verify that this coincides, in the limit k\ <C k 2 ,k%, with the bispectrum for scalar field perturbations given 
by Seery and Lidsey in Ref. I16j for the case of a scale-invariant spectrum. 



C. Non-Gaussianity of curvature perturbations 

The power spectrum of the curvature perturbation due to inflaton fluctuations during single field inflation is given 
by the first-order relation $Zffy . Substituting this into Eq. (|T%|l gives 

H ,2 



p d*)=\jr) P Q( k )- (48) 

The power spectrum of scalar field perturbations at Hubble exit is given by Eq. (|37|) and thereafter the non-linear 
curvature perturbation remains constant, on super-Hubble scales, for adiabatic perturbations even after slow-roll 
inflation ends 29]. Thus for super-Hubble modes with k <c aH we can write the power spectrum in terms of 
quantities at Hubble-exit 

7> c (*) s ^P c (fc)=(|L) . (49) 
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Slow-roll gives an almost scale-invariant power spectrum of curvature perturbations on super-Hubble scales, whose 
tilt is given by 

n-1 = ^i = -6e + 2 V . (50) 
a Ink 

We can write the bispectrum 1|17|) for the curvature perturbation in terms of a local transformation between second- 
order gauge- invariant variables, Eq. I|31|l . and the intrinsic bispectrum of the scalar field perturbations. Substituting 
Eqs. (J2HJ) and into Eq. JTHJ) we obtain 

3 



B c (k 1 ,k 2 ) = -(^-j B Q (k 1 ,k 2 ) + (2e-r ] )(P ( (k 1 )P c (k 2 )+P < (k 2 )P < (k 3 ) + P c (k 3 )P c (k 1 )) . (51) 

We see that if the field perturbations on spatially flat hypersurfaces were strictly Gaussian at second order, Bq = 0, 
then the second order curvature perturbation ( 2 would be of the "local" form defined in Eq. (121 H with /nl — 5(r;— 2e)/6. 
In the squeezed triangle limit, k <C the scalar field bispectrum at Hubble exit is given by l|47|) and thus 

B c (k, k,) = (6e - 2 V )P ( (k)P c (K) , (52) 

where we have used P((k) 3> -Pf(fc*) for a nearly scale-invariant spectrum. Thus comparing with Eq. (|21|) in this 
squeezed triangle limit, we can identify 

f NL = ^(-6e + 2r 1 ). (53) 

As noted by several authors 0,0] this gives a consistency relation between the primordial non-Gaussianity and the 
tilt of the power spectrum l|50|) 

fNL = ^(n-l). (54) 

which does not rely on slow-roll during inflation, but does rely on the adiabaticity of the perturbations on super-Hubble 
scales and thus the existence of a conserved curvature perturbation after Hubble-exit. 



IV. MULTI-FIELD INFLATION 



This discussion of of the non-Gaussianty of primordial curvature perturbations from inflation can be readily extended 
to the multi-field [22] case by generalising Eq. i|32|) 

C = N ^ + N ,iJ s Qi s Qj + ■■■ ■ ( 55 ) 

i i,j 

where we identify £ with the perturbed expansion of a uniform density hypersurface with respect to an initial uniform- 
expansion hypersurface during inflation, and we define = dN/difi. Hence to first and second order we have 

i 

(2 = J2 N ' iS2Qi+1 12 N >v SiQiSiQ i ■ ( 5? ) 

i i,j 

The gauge- invariant definition of field perturbations on uniform-expansion hypersurfaces, 5Q, was given in Eqs. I|24l) 
and 

It is important to realise in the multi-field case that the primordial curvature perturbation is evaluated on a uniform- 
density hypersurface at the some primordial epoch after inflation has ended. Thus the integrated expansion to this 
uniform-density hypersurface as a function of the initial field local values, N((pi), must in general be calculated not 
only during inflation, but throughout the subsequent cosmological history including reheating. Only in the case of 
adiabatic perturbations on super-Hubble scales, where the perturbed evolution follows the background trajectory in 
phase space |3lll33| . can we equate SN with the perturbed expansion at Hubble-exit, independently of the subsequent 
expansion history, as was done in Eqs. I|32H3()|I for the single-field case. Nonetheless, with this caveat we can formally 
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write down the primordial perturbation in terms of the integrated expansion N and it's dependence on the different 
fields. 

Thus substituting Eqs. I|56(l and (|57Jl into the definitions of the power spectrum l|18ll and bispectrum (jf 9(1 . we obtain 
P c (k) = ^NjNjGjik), (58) 

3 

B c (k u k 2 ,k 3 ) = ^ A., A., A, ,,/.';,' '''.:/,,./-'•../,•;) • N ,i N J N ,™ [C im (ki)C jn (k 2 ) + perms] , (59) 

where we define the multi-field cross-correlation and bispectrum at leading order 

(6iQi*8iQ jy ) = Cy(fc) ^ 3 (k + k'), (60) 
{SQ i>kl SQ jtka SQmMa) = B%- m (ki, k 2 ,k 3 ) ^ 3 (k x + k 2 + k 3 ) . (61) 

Equation l|59|) allows for an arbitrary bispectrum, Bq 1 '" 1 , for the field perturbations on spatially flat hypersurfaces 
and generalises the result of Lyth and Rodriguez 22] to allow for correlations between the fields on super-Hubble 
scales. In principle it does not rely on slow-roll, only that the subsequent expansion can be expressed as a function 
solely of the local field values. Nonetheless in the rest of this paper we will restrict our analysis to slow-roll in order 
to calculate the bispectrum of the field perturbations. 



A. Slow-roll field perturbations 

The field perturbations originate as vacuum fluctuations on sub-Hubble scales and hence SQi and SQj are indepen- 
dent random variables at Hubble-exit with 

jj 

to leading order in a slow-roll expansion, where 

8 

47rfc 3 



ei,kSi,k'> = -rb S^(k + k'). (63) 



At Hubble-exit we have Cy(A;*) = (H 2 /2k^)8ij, but coupled evolution of linear perturbations on large scales (propor- 
tional to d 2 V/dtpidfj in slow-roll approximation) can give rise to off-diagonal terms in multi-field inflation plL |42|. 
Such off-diagonal terms are small in the slow-roll approximation, but their effect can become significant when inte- 
grated over many expansion times for long- wavelengths, <C 

Once again we can estimate the bispectrum of the scalar field perturbations in the squeezed triangle limit k -c fc* 
at Hubble-exit k* = aH, due solely to gravitational coupling in slow-roll inflation. Generalising Eq. (|42|l to multi-field 
inflation we obtain 

^8 2 Qi,k* = ^2~^-{8iQj •8 1 Q % ) k „ , (64) 



H 

j 



where Hi = dH/difi. Hence for k\ <C k 2 — k 3 ~ k* 



(<5Qi,ki 8Qj t k 2 SQm.ks) = X! "cT ({fiiQiM (8iQn • 8iQj)k 2 £iQ m ,k 3 } + (SiQiM #iQj,k a (8iQ n • £iQm)k 3 )) , 

n 

= 2 ^T C in(h)C jm (k 2 ) # 3 (k ± + k 2 + ks) . (65) 

n 

Comparing l|t)5[) with the definition of the multi-field bispectrum in Eq. l|61f) we see that in the squeezed triangle 
limit k <C fc* at Hubble-exit 

B^ m (k,K) = Y^C m (k)C 3m (K). (66) 



As we would expect, the bispectrum for field fluctuations due to gravitational coupling is non-zero only when the 
Hubble rate at Hubble-exit is dependent upon the field values. 
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Following Ref. [33 we can use the unperturbed background solution, <fii(t), to define the instantaneous adiabatic 
or "inflaton" perturbation lying along the background trajectory in multi-field inflation. 

da = — <j9 ' dtfj . (67) 

All orthogonal directions in field-space describe instantaneous isocurvature or entropy field perturbations 0, 03 which 
we denote generically by 5x- 

In the slow-roll approximation the background trajectory is given by the gradient of the potential, or equivalently 
the Hubble rate. Thus 

£ ^iQn = ^fhQa ■ (68) 

n 

Thus the bispectrum (|66|) is only non-zero for long-wavelength field perturbations SQi which are correlated with the 
instantaneous inflaton direction, a: 

/>';;/"";/■••/.•.; = -2e^-C itr {k)C jm (K) , (69) 

analogous to the single- field case (|47|) . where the slow- roll parameter (|29[l is given by 

dH & 

Isocurvature field perturbations that remain decoupled from, and hence uncorrelated with, the inflaton perturbation 
(i.e., C' x<T (k) — 0), retain a purely Gaussian distribution with B x (k,k*) = unless we consider non-gravitational 
interactions, neglected here in our slow-roll approximation. 



B. Non-Gaussianity of curvature perturbation in slow-roll 

Finally we can write the bispectrum l|59|l for the curvature perturbation for multi-field slow-roll inflation in the 
squeezed triangle limit, k <C k*, as 



B ( (k, K) = 2 NiN tj ( eN d N a C itT {k) + £ N jn C in (k) j P Qj (K) . (71) 

One verify that this reduces to Eq. (|52|l in the single field case where we have 7V CT and N^ aa given by Eqs. 135fl and 
(J2HJ). This single- field result for the non-Gaussianity still holds in a multi-field setting so long as the expansion history 
is independent of all the other fields. 

An alternative limit is the case where the expansion history is dominated by an isocurvature field, such that 
|iV x | 3> \N t &\. This would apply for instance to the curvaton scenario |43|. In this case the primordial curvature 
perturbation is given by 

P c (k) = (N x fP x (k). (72) 

There is no non-Gaussianity due to gravitational coupling, B x = 0, for an isocurvature field that remains decoupled 
from the from the inflaton, so that C xa — 0. Hence the non-Gaussianity of the primordial curvature perturbation 
(j59*|l is given by 

B ( (ki,k 2 ,k 3 ) = (N X ) 2 N XX [P x (fc x )P x (fc 2 ) + perms] , 
N 

= [P C (ki)P C (k2) + perms] . (73) 

This is purely of the local form Q21H and we can identify [2^] 

6 (N a y 

Even in the general case it can be shown that the contribution of the intrinsic non-Gaussianity of the field pertur- 
bations to the bispectrum of the curvature perturbations must be small in the slow- roll approximation |24j . 
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V. CONCLUSIONS 



We have shown how the non-Gaussianity of primordial curvature perturbations is described using gauge-invariant 
cosmological perturbations at second-order. The split between zeroth, first and second-order parts of inhomogeneous 
fields is arbitrary, however if we expand about a FRW cosmology, the zeroth-order part of cosmological fields can be 
taken to be spatially homogeneous. In this paper we have further assumed that the first-order part can be taken to be 
a Gaussian random field, and the second-order part describes the non-linear evolution of inhomogeneities, and hence 
the non-Gaussian part of the distribution. 

The arbitrariness of the choice of coordinates leaves a residual gauge-freedom. Nonetheless one can construct 
gauge-invariant quantities, at any order, by specifying an unambiguous, physical quantity [3^ |. In particular we have 
given a manifestly gauge-invariant definition of £, the perturbed expansion on uniform density hypersurfaces on large 
scales, used by many authors to describe the primordial metric perturbation as it remains constant, at any order, for 
adiabatic perturbations [2^, I30L I35L l36| . This is related to the dimensionless density perturbation on spatially flat 
hypersurfaces, however the non-linear relation between the two leads to a non-Gaussian metric perturbation for a 
purely Gaussian density perturbation, and vice versa. 

In the case of single-field inflation, the field perturbations become adiabatic in the large-scale limit 01 an d the 
primordial C long after inflation has ended can be equated with the perturbed expansion on uniform-field hypersurfaces 
at Hubble exit. We have shown how this is equivalent to the non-linear extension of the S N- formalism [22l I2H l28l 
l3(l | which allows one to calculate the primordial perturbation in terms of SQ, the field perturbations on spatially 
flat hypersurfaces, at Hubble exit and the unperturbed FRW solution for the integrated expansion, N((p). This 
formalism is especially compact for multiple fields where the primordial metric perturbation (say, during primordial 
nucleosynthesis) may be completely different from that at Hubble exit during inflation. On the other hand one should 
bear in mind that N(tpi) is a non-linear function of several scalar fields which must be determined by solving the 
background field equations to determine the integrated expansion up to some fixed primordial density as a function of 
the field values during inflation. Nonetheless Lyth and Rodriguez |22( have recently shown how powerful this method 
is for describing the non-Gaussianity of the primordial perturbation due to non-linear evolution on super-Hubble 
scales. 

In single- or multi-field inflation one can give the leading-order expression for the bispectrum of the primordial 
metric perturbation in terms of the power spectrum and the intrinsic bispectrum of the fields at Hubble-exit during 
inflation. We have used a simple argument to estimate the intrinsic bispectrum of the field perturbations at Hubble 
exit due to gravity in the limit of degenerate triangles during slow-roll inflation with single or multiple fields, and hence 
give an expression for the non-linearity of the primordial density perturbations from slow-roll inflation in this limit. 
In line with previous results we find that the non-Gaussianity of field perturbations is suppressed during slow-roll, 
and we have shown that it vanishes at leading order for isocurvature perturbations during inflation. Our final result 
in Eq. (|71|l generalises earlier work to allow for non-scale invariance of the scalar field power spectra and correlations 
between the fields that can develop on super-Hubble scales. 
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